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ABSTRACT 



ISO observations at 200/im have modified our view of 
the dust component in spiral galaxies. For a sample 
of seven resolved spirals we have retrieved a mean 
temperature of 20K, about 10K lower than previ- 
ous estimates based on IRAS data at shorter wave- 
lengths. Because of the steep dependence of far- 
infrared emission on the dust temperature, the dust 
masses inferred from ISO fluxes are a factor of 10 
higher than those derived from IRAS data only, lead- 
ing to gas-to-dust ratios close to the value observed 
in the Galaxy. The scale-length of the 200/im emis- 
sion is larger than for the IRAS 100/im emission, with 
colder dust at larger distances from the galactic cen- 
tre , as expected if the interstellar radiation field is the 



ma in source of dust heating. The 200/im scale- length 
is also larger than the optical, tor all the galaxies m 
the sample. This suggests that the dust distribution 
is more extended than that of the stars. 

A model of the dust heating is needed to derive the 
parameters of the dust distribution from the Far- 
Infrared (FIR) emission. Therefore, we have adapted 
an existing radiative transfer code to deal with dust 
emission. Simulated maps of the temperature distri- 
bution within the dust disk and of the dust emission 
at any wavelength can be produced. The stellar spec- 
tral energy distribution is derived from observations 
in the ultraviolet, optical and near infrared. The pa- 
rameters of the dust distribution (scale-lengths and 
optical depth) are chosen to reproduce the observed 
characteristics of the FIR emission, i.e. the shape of 
the spectrum, the flux and the spatial distribution. 
We describe the application of the model to one of 
the galaxies in the sample, NGC 6946. 



is hotter and emits preferentially at shorter wave- 
lengths. The link between extinction and cold dust 
is revealed by the similar values of the gas-to-dust 
mass ratio retrieved in the Galaxy from independent 
studies of e xtinction and emission. Using mean dust 
properties (Hildebrand 1983;Whittetl992), a value of 
130 can be derived Irom the correlation between the 
hydrogen column d ensity and the B-V colour excess 
( [Bohlin ct al. 1978 ). From the Galactic FIR emission 
at A > 100/im observed b y the instrument DIR BE 
aboard the satellite COBE, Sodroski et al. (1994) de- 
rived a mean value of 160. 

Until recently, the main source of FIR data for spiral 
galaxies has come from the IRAS satellite. Because 
of its limited spectral coverage ( A < 120/iin). IRAS 
was no t able to detect cold dust. Devereux fc Young 



1990) measured the gas-to-dust ratio for a sample of 
8 spiral galaxies, using the IRAS 60 /im and 100/im 
fluxes. A mean value of 1080 was derived, almost an 
order of magnitude larger than in the Galaxy. As- 
suming that the Galaxy is an average spiral, they ex- 
plained the discrepancy with 90 per cent of the dust 
mass being at T ~ 15 K, too cold to be detected by 
IRAS. 

Therefore, observations at longer wavelengths are es- 
sential to detect cold dust and to trace the internal 
extinction in a spi ral galaxy. In this paper we de- 
scribe ISOPHOT ( |Lemke et al. 19961 ) observations 
at 200 /tm of a sample of seven spiral galaxies. An 
extended cold dust emission has been revealed (Sec- 
tion EJ) . A model of FIR emission has been built to 
derive the parameters of the dust distribution from 
the observed properties of dust emission (Section pi). 
A summary and a discussion of the results are given 
in Section |4~| 
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2. ISO OBSERVATIONS OF SPIRALS 



1. INTRODUCTION 



Diffuse dust is responsible for the internal extinction 
in a spiral galaxy. The grains of this diffuse com- 
ponent are heated up to temperatures T < 20 K by 
the mean galactic Interstellar Radiation Field (ISRF; 
Reach ct al. 1995J ). Emission from cold dust peaks 
at A > 150/iin, while dust in circumstellar regions 



We give in this section a summary of the ISOPHOT 
200 /im observ ations and results. A full description 
is presented in Alton ct al. (1998)|. 



Eight objects were observed with ISOPHOT, seven 
spiral galaxies and one irregular. The targets were se- 
lected to be representative of quiescent objects with 
normal FIR-to-blu c luminosity ratio (Lfj.r/ Lb < 1, 
3oifer et al. 1987) and to have large apparent sizes 
(~ 10 arcmin), to be resolved by the ISO beam 
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Figure 1. Images of the spiral galaxy NGC6946, 
Alton et al. (1998) 



in the B-band (left) and at 200 [im (right). Images are from 



(FWHM=117 arcsec, [Tuffs ct al. 1996[ ). B-band 
images were also obtained, together with resolution 
enhanced HiRes IRAS images at 60 and 100 /im 



QR|ccl993|). The ISOPHOT images are similar in 
morphology to the lUU/im HiKes, that have a com- 
parable resolution. IRAS and B-band images were 
convolved to the ISO resolution and registered with 
the 20Qjim images. As an example, we present in 
Figure l] a full-resolution B-band image and a 200 /im 
image for one of the objects of the sample, the spiral 
galaxy NGC 6946. 

Total integrated fluxes at 60, 100 and 200 /im were 
used to derive dust temperatures and masses. As for 
the dust emissivity, the form Q em (A) = Qem(Ao) x 
(A/Ao) - ^ was chosen, allowing for two different wave- 
lengt h dependences. 8= \ and 2, and using Q em (Ao) 



from Hildcbrand (1983). After computing the dust 



mass, gas-to-dust mass ratios were derived usin g the 



tot al (atomic+molccul ar) gas masses listed by Dev 



ere |ux & Young (1990) 
carried out 



the IRAS 



Two set of calculation were 
view, using only the 60 and 
100 /im fluxes from IRAS; and the ISO view, using 
the 100 /xm IRAS flux and the 200 /im ISO flux. The 
mean values of dust temperatures and gas-to-dust 
mass ratios over the observed sample are presented 
in Table g. 

Clearly, a larger amount of cold dust is revealed using 
the 200 /im data. Dust temperatures from the ISO 
view are colder by almost 10 degrees, and the gas- 
to-dust ratios are smaller by one order of magnitude, 
with respect to the IRAS view. The values for the 
gas-to-dust ratio are similar to the Galactic. The cali- 
bration of ISOPHOT is still uncertain. An integrated 
flux at 200 /im was derived for N GC 6946 from the 
observations of Engargiola (1991). The ISO flux for 
the same galaxy is larger by 30 per cent. The same re- 



sult is obtained comparing the observed 200 /im back- 
ground with that extrapolated from 100 urn, using a 



Galact ic spectrum at high latitude from Reach et al. 



(1995). However, even accounting for an over estima- 



tion Oi the ISOPHOT calibration (results in brackets 
in Table nj) , a larger quantity of cold dust is revealed 
by the ISO view. 

Uncertainties on the calibration do not affect the re- 
sults about the spatial distribution of FIR emission. 
Azimuthally averaged profiles were obtained for the 
optical and FIR smoothed images, and the exponen- 
tial scale-lengths were measured between 1.5 arcmin 
and 3.5 arcmin. The mean scale-lengths for the seven 
spirals in the sample are presented in Table 0, nor- 
malised to the 200 /jm value. The larger 200 /im scale- 
length with respect to the other FIR images shows 
that dust is colder at larger radii, consistently with 
dust grains heated by a diffuse ISRF of higher in- 
tensity in the centre of the galaxy. The 60 /im scale- 
length is similar to the 100 /im, a result also observed 



Table 1. Mean temperatures and gas-to-dust ratios 
for the observed sample, derived for the IRAS view 
and the ISO view, for two choices of the emissiv- 
ity spectral index (3. Values in brackets allow for an 
over estimation of the ISO fluxes by 30 per cent. 



P IRAS 



ISO 



(T(K)> 
(Gas-to-Dust) 



34 21 (24) 

28 18 (19) 

2800 230 (494) 

5100 220 (446) 



Table 2. Mean scale-lengths of the sample 



60/200 0.57±0.05 
100/200 0.56±0.04 
B/200 0.79±0.05 



in the Galaxy flSodroski ct al. 198E ) and explained 
with emission from small grams whose heating con- 
ditions do not depend heavily on the ISRF gradient. 
T he most striking results is the smaller B band scale- 



length with respect to that at 200 /im. If the dust 
distribution is similar to the stellar, the scale-length 
of FIR emission should be smaller than the optical, 
because the dust temperature is higher in the cen- 
tre of the disk. Therefore, the larger 200 /mi scale- 
length could be a hint for a dust distribution more 
extended than the stellar. Extended dust distribu- 
tions also emerge from models of COBE e mission at 
140 /im and at 240 /im ( |Davics ct al. 1997 ) and from 
fit s of surfa ce brightness in edge-on galaxies ( Xilouris 
et kl. 1999|). 



To check tor the influence ot transient ettects along 
the scan direction, the galaxies were divided in four 
quadrants and the scale-lengths analysis repeated 
only on the two quadrants perpendicular to the scan 
direction. The same result was obtained. 



3. MODELLING THE FIR EMISSION 



In this section we briefly describe a model of FIR 
emission in spiral galaxies. The main purpose of the 
model is the derivation of the geometrical and optical 
parameters of the dust distributions through a com- 



parison oi simulations with observed. FIR emission. 



In pailiuulai, wu wanled Lu Lest Ihe hypulliusis lliaX 
the extended dust emission described in Section Ej is 
due to an extended distribution of dust. A detailed 
description of the mo del and its results is deferred t o 
a forthcoming paper ( Bianchi, Davies fc Alton 1999| ). 



3.1. Description of the model 



The model is based on a simplified v ersion of t,ho i 
Monte Garlo radiative tra nsfer code of Bianchi, Fer 



rara fc Giovanardi (1996) . The Monte Carlo method 
allows for a correct treatment of multiple scattering 
in geometries typical of spiral galaxies. For the mod- 
els presented here, we use a double exponential disk 
geometry for both stars and dust, although the code 
is able to work with other distributions. The stellar 
and dust disk are geometrically defined by the radial 
and vertical scale-lengths. The opacity of the dust 
disk is defined by the face-on optical depth through 
the centre, Ty. We used em pirical dust properties for 



absor ption and scattering (|Gordon, Calzetti fc Witt 



1997) 



pir 



The code is monochromatic. Simulated images for a 
particular wavelength are produced, together with a 
map of the energy absorbed by dust from starlight at 
that wavelength. Therefore, to obtain a map of the 



total energy absorbed by dust from the whole stel- 
lar Spectral Energy Distribution (SED), simulations 
are produced for 17 wavelengths bands covering the 
spectral range of stellar emission, from the ionization 
limit to the Near-Infrared. 

The main aim of this work is the modelling of the FIR 
emission, that is produced by grains heated at ther- 
mal equilibrium. Since part of the absorbed energy is 
re-emitted in the Mid- Infrared by transiently heated 
small grains, a correction is applied to subtract their 
contribution from the absorbed energy map. The 
correction is derived from the abso rption efficiency of 
small grains in t he dust model of Desert, Boulanger 



fc Puget (1990) 



From the corrected absorbed energy map the tem- 
perature distribution along a meridian plane of the 
galaxy model is retrieved. A key quantity for this 
calculation is th e cmissivitv. We have derive d the 



emissivity value (Bianchi, Davies & Alton 1999) from 
maps of Galactic FIR emission, temperature and ex- 
tinction (Schlegel, Finkbciner & Davis 1998), assum- 
ing a spectral dependence as measured on high signal- 
to-noise FI R spectra of the Galactic plane (Reach 



et al. 1995) 



Images of FIR emission are then produced from the 
temperature maps and dust distribution, integrating 
along specific line of sights. In this paper we describe 
an application to the spiral galaxy NGC 6946. Sim- 
ulated optical and FIR images have been produced 
f or the inclination of 34° . appropriate for the g alaxy 



(Athanassoula, Garcia-Gomez & Bosma 1993), for 
several star-dust config urations. W e have adopted 
a distance of 5.5 Mpc flTully 1988 ) for the galaxy; 
an intrinsic unextinguished stellar radial scale-length 
a* = 2.5 kpc (derived from a K-band image); a 
radial/vertical scale-lengths ratio a*//3* = 14.4 (a 
mean value derived from the compilation of vertical 
scale-lengths for different stellar types of Wainscoat 



et al. (1992)|) 



3.2. Results 

We first tried a model with a standard distribution 
for the dust disk, i.e. with the dust radial scale- 
length of the same dimension as the stellar (a<j = a 4 ) 
and vertical scale-length one half of the stellar ((3d = 
0.5/3*). The choice of the dust vertical scale-length 
is mainly dictated by the observation of extinction 
lanes in edge-on galaxies, that cannot be explained 
with dust distribution thicker than the stellar. 

The SED of the stellar and dust emission are pre- 
sented in Figure H for standard models of different Ty. 
The optical data-points for the stellar emission are 
from Engargiola (1991) and represent the flux inside 
the B-band half-light aperture (5 arcmin diameter). 
UV fluxes inside the same aperture a re derived from 
Rifatto, Longo fc Capaccioli (1995) . FIR fluxes in- 
side the same aperture are derived from IRAS HiRcs 
images an dfrom the 200 /im I SO image. The 160 /im 
flux from Engargiola (1991)| has been included as 
well. The simulated optical and UV images were in- 
tegrated inside the B-band half light radius (also de- 
rived from the simulation) and the fluxes normalised 
to the observed values. The solid line in Figure @ 
represent this SED of the observed stellar emission. 
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Figure 2. Optical and FIR SEDs for standard models 
(ad = a*; $d = 0.5/3*/ of different optical depths. 



The SED of the FIR emission (and that of the intrin- 
sic, unextinguished, stellar radiation) are presented 
for models of different optical depths. 

The SEDs of FIR radiation in Figure show that 
only optically thick models (ry ~ 5) can provide the 
observed amount of energy absorbed (and emitted) 
by dust. As for the scale-length of FIR emission, 
they increase with the optical depth. However, the 
B hand scale-length increases with extinction as well. 



Be pause of this, the smaller value for the B/200 scale- 
lengths ratio occurs for optically thin model with 
= 0.5 (1.4 vs 0.9 measured for NGC 6946V The 



TV 

ter iperature distributions arc quite similar for any of 
the dust disk models. The temperature distribution 
on the meridian plane for the standard model with 
Ty = 5.0 is presented in the top panel of Figure Q 
The model temperature s are compatible wi th those 
observed in the G alaxy (jSodroski et al. 1997J ) and on 
NGC 6946 itself (|D^ie^ 

I I 

Xilouris et al. (1999) modelled the surface brightness 
of a sample of seven edge-on galaxies with a double 
exponential geometry for dust and stars. The ex- 
tinction lane is well fitted if the dust distribution is 
thinner and has a larger radial scale-length than the 
stellar ((ad/a*) = 1.5, ((3 d/ (3*) = 0.5). The dust disk 
are found to be optically thin, with a mean Ty = 0.5. 
Figure [| shows the result for the case with extended 
disk. Despite the internal extinction has increased, 
still optically thick models are needed to absorb the 
right amount of energy. In an extended model, larger 
amounts of colder dust are present in the outer skirts 
of the disk (Figure 5|) and this causes an increase in 
the FIR scale-lengths. However, the increase in the 
B-band scale-length due to extinction compensates 
for the larger FIR scale-lengths and the B/200 scale- 
lengths ratios are larger then observed. The smaller 
value (1.15) is again achieved for the optically thin 
case with Ty = 0.5. The models of COBE emission of 



Figure 3. Optical and FIR SEDs for extended models 
(ad = 1.5a*, (3d = 0.5/3*/ of optical depths Ty = 0.5 
and Tv = 5.0. The standard model with tv — 5.0 is 
also shown, for comparison. 



Davies et al. (1997)| suggest that the dust distribution 
is more extended then the stellar also in the vertical 
direction ((3d = 2/3*). The inclusion of a thicker dust 
disk, however, does not modify the trend seen for the 
other models. 

Since the molecular gas distribution of NGC_6JH£ 
follow closely the emission detected by IRAS, ~~ 



Davies 



et al. (1999)| suggested that the broader 200 /im pro 



file could be due to a more extended dust distri 
bution associ ated with the atomic gas. 



Tacconi fc 
Young (1986)| report the column density or molecular 
and atomic hydrogen as a function of the galactocen- 
tric distance. The molecular gas has a steep profile, 
with an exponential scale- length wla 4 . The atomic 
gas distribution has a dip in the centre, reaching 
a maximum at «2a*, then declining exponentially 
with a scale- length «3a*. The atomic and molecu- 
lar phases have the same mass, but the former has a 
much broader distribution. A model with two dust 
disk mimicking the gas distributions is presented in 
Figure 0. As central face-on optical depth we choose 
Ty—5 for the dust associated with the molecular gas, 
and t v =5 for the dust associated with the atomic 
phase. The adopted optical depths were derived us- 
ing the correl ation between ext inction and gas col- 
umn density ( Bohlin et al. 1978) and a st andard ex- 
tinction law for the Galaxy ( Whittetl992 ) . 



Despite the large extent of the disk associated with 
the atomic gas, the FIR emission in the model is 
dominated by the optically thick disk associated with 
the molecular phase. The behaviour of the double 
disk model is therefore similar to that of a standard 
model with tv—5. A good fit can be found only for 
the SED while the B/200 scale-length ratio is larger 
than observed. Interestingly, a better match of both 
FIR emission and scale-lengths ratio can be found 
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Figure J h Optical and FIR SEDs for the double disk 
model, and for a single disk model with ad — 3a* and 
Tv =5.0. The standard model with ry = 5.0 is also 
shown, for comparison. 



for a single disk model, of the same scale-length of 
the one previously used for the atomic gas. An op- 
tical depth ry =3 is needed to match the observed 
fluxes (Figure H). The bottom panel of Figure [~ 
shows the temperature distribution for this mode 
increasing the dust scale-length with respect to the 
stellar results in larger amounts of colder dust in re- 
gions of low intensity ISRF. The FIR scale- length has 
increased enough to compete with the change in the 
B-band scale-length due to extinction and the result- 
ing B/200 scale-length ratio is now slightly smaller 
than 1. 



4. SUMMARY & CONCLUSIONS 



ISO observations of a sample of seven spiral galax- 
ies at 200 fjm have revea led large amounts of cold 
dust ( [Alton et al. 1998| ). The FIR emission is 
found to have a larger scale-length than the optical 
({otB / &200) ~ 0-8), thus suggesting a broader distri- 
bution for dust than for stars. 



We produced a model for the radiative transfer 
and FIR emission for geometries typical of a spiral 
galaxy. We applied the model to the observations of 
NGC 6946. For all the model explored, an optically 
thick dust distribution (with central face-on optical 
depth Ty ~5) is necessary to match the FIR spec- 
trum. However, optically thick models fail to repro- 
duce the observed ratio of B-band and 200 ;Uiri scale- 
lengths. For models with dust radial scale-length 
larger than the stellar by a factor 1.5, as suggested 
by fitting of surface brightness in edge-on galaxies, 
the optical-FIR scale-lengths ratio is closer to the 
observed value only in optically thin cases. A better 
fit of both FIR spectrum and spatial distribution of 



emission, in the optically thick case, can be reached 
only for O^/a* > 3. 

A few approximations have been adopted in this 
work. We have tested that the results are non qual- 
itatively modified by: the neglect of the ionising 
stellar continuum; the presence of a small stellar 
bulge; the assumption of a single stellar vertical scale- 
length, instead of having a thinner stellar d isk for 
young objects ( Bianchi, Davies fc Alton 1999J ) emit- 
ting preferentially at smaller A. We have also as- 
sumed a single stellar radial scale-length at any A, 
following the idea that colour gradients are mainl y 
produced by du st extinction ( Peletier et al. 1995|) . 
De Jong (1996) claims that extinction cannot explain 
the colour gradients observed in a sample of 86 spiral 
galaxy, for reasonable models of galactic disks. How- 
ever, this definitions does not seem to include models 
with dust scale-lengths larger than the stellar, that 
have colour gradients as large as the observed. 

A major assumption in the radiative transfer model 
is that of a smooth distribution, both for stars and 
dust. Clumping may be able to explain the paradox 
of needing both optically thick and optically thin dis- 
tributions, to reproduce the FIR spectrum and scale- 
lengths ratio, respectively. Dust associated with the 
gas in molecular clouds may be responsible for most 
of the FIR emission, if large amounts of energy are 
absorbed from embedded stars. On the other hand, 
an optically thin distribution correlated to the atomic 
gas will be able to explain the observed spatial dis- 
tribution of FIR emission. The galaxy may look as 
optically thin, when seen face-on. as derived fro m the 
surface brightness fits of Xilouris et al. (1999) 



Such a scenario, however, does not seem to be sup- 
ported by observations. For large objects with a 
well resolved clump structure, like the Milky Way 
and M31, only a 30 per cent of the FIR emission 
comes from star- forming regions embedded in molec- 
ular clouds ( ISodroski et al. 1997J ; [Xu fe Helou 1996| ). 
It is then difficult to evaluate the effect ol clumping 
on extinction and FIR emission. A study of the in- 
fluence of clumping on the radi ative transfer through 
a galactic disk is presented by Bianchi et al. (1999). 
It is found that clumping affects only marginally the 
radiative transfer. It would therefore be difficult for 
clumpy dust associated with molecular clouds to pass 
"undetected" in surface brightness fits. However, re- 
sults depend highly on the description adopted for 
the clumps distribution. A description appropriate 
to the molecular distribution of NGC 6946. rather 
than the Galactic adopted in Bianchi et al. (1999)1 , 
may produce different results. A proper model of FIR 
emission for a clumpy dust distribution is therefore 
needed to confirm the prediction of the last para- 
graph. 
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